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^ ■ ABSTRACT 

C*~) ' Using a series of high-resolution N-body simulations of the concordance cosmology we 

04 , investigate how the formation histories, shapes and angular momenta of dark-matter 

haloes depend on environment. We first present a classification scheme that allows to 
distinguish between haloes in clusters, filaments, sheets and voids in the large-scale 
distribution of matter. This method (which goes beyond a simple measure of the local 
density) is based on a local-stability criterion for the orbits of test particles and closely 
\ relates to the Zel'dovich approximation. Applying this scheme to our simulations we 

■ then find that: i) Mass assembly histories and formation redshifts strongly depend on 

— ' environment for haloes of mass M < M* (haloes of a given mass tend to be older in 

\Q , clusters and younger in voids) and are independent of it for larger masses (M* here 

' indicates the typical mass scale which is entering the non-linear regime of perturbation 

growth); ii) Low- mass haloes in clusters are generally less spherical and more oblate 
than in other regions; Hi) Low-mass haloes in clusters have a higher median spin than 
I 1 in filaments and present a more prominent fraction of rapidly spinning objects; we 

identify recent major mergers as a likely source of this effect. For all these relations, 
+3 . we provide accurate functional fits as a function of halo mass and environment. We 

£2 ' also look for correlations between halo-spin directions and the large-scale structures: 

. . | the strongest effect is seen in sheets where halo spins tend to lie within the plane of 

symmetry of the mass distribution. Finally, we measure the spatial auto-correlation of 
k> ' spin directions and the cross-correlation between the directions of intrinsic and orbital 

angular momenta of neighbouring haloes. While the first quantity is always very small, 
' we find that spin-orbit correlations are rather strong especially for low-mass haloes in 

clusters and high- mass haloes in filaments. 
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1 INTRODUCTION 

Numerical simulations and analytical work have shown that 
the gravitational amplification of small density fluctuations 
leads to a wealth of structures resembling the observed large- 
scale distribution of galaxies. The resulting mass density dis- 
tribution can be thought of as a "cosmic web" l|Bond et al.l 
1996) characterised by the presence of structures with differ- 
ent dimensionality. Most of the volume resides in low-density 
regions (voids) which are surrounded by thin denser sheets of 
matter. A network of filaments of different sizes and density 
contrasts departs from the sheets and visually dominates 
the mass distribution. Dense clumps of matter lie at the in- 
tersections of filaments. From the dynamical point of view, 
matter tends to flow out of the voids, transit through the 



sheets and finally accrete onto the largest clumps through 
the filaments. 

In a Universe dominated by cold dark matter (CDM), 
this description applies only after coarse-graining the density 
distribution on scales of a few Mpc. On smaller scales, the 
power in the primordial spectrum ends up producing a hier- 
archical distribution of (virialised) dark-matter haloes whose 
positions trace the large-scale structure described above. Ac- 
cording to the current cosmological paradigm, galaxies form 
within these haloes. 

Astronomical observations show that galaxy properties 
in th e local Univers e vary systematically wi t h environment 
(e.g. iDressleJ ll98Cll : [Kauffma nn et al. 1 12004 iBlanton et al.l 
2005). As a fundamental step towards understanding galaxy 
formation it is thus important to establish how the prop- 
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erties of dark-matter haloes depend on the environment in 
wh ich they reside. A first a t tempt in this direction was made 
by iLemson fc Kauflhiannl (|l999T j who found that mass is 
the only halo property that correlates with environment 
at variance with concentration, s pin, shape and forma - 
tion epoch. Using marked statistics, ISheth fc Tormenl (|2004l ) 
found evidence that haloes of a given mass form earlier 
in dense regions. Higher resolution simulations confirmed 
this finding and helped to better quantify it as a function 
of halo mass and redshift (iGao et alj 20051; Croton et al 



20061 ; lHarker et alll2006l ; iReed et al.ll2006l ; lMaulbetsch et al 
2006). At the same time it has become clear that also other 
halo properties as concentration and spin correlate with 



local environmen t llAvila-Reese et al 



2005 1; iBett et al.l 120061 ; iMaccio et al 
2006n . 



2005; 



Wcchsler et al 



2006: Wetzel et al 



Although the large-scale structure of matter is promi- 
nently reflected in the halo distribution, no efficient auto- 
mated method has been proposed to associate a given halo 
to the dynamical structure it belongs to. Most of the en- 
vironmental studies mentioned above use the local mass 
density within a few Mpc as a proxy for environment. In 
this paper we follow a novel approach and associate dark- 
matter haloes to structures with different dynamics. Voids, 
sheets, filaments and clusters are distinguished based on a 
stability criterion for the orbit of test par ticles which is in - 
spired by the Zel'dovich approximation (|Zerdovichlll97(]| ). 
Our method is accurate, fast, efficient and contains only one 
free parameter which fixes the spatial resolution with which 
the density field has to be smoothed (as in the evaluation of 
the density). We show that any classification based on local 
density is degenerate with respect to ours which we regard 
as more fundamental. We find that all halo properties at 
zero redshift show some dependence on the dynamical envi- 
ronment in which they reside. We accurately quantify this 
dependence and show that halo properties smoothly change 
when one moves from voids to sheets, then to filaments and 
finally to clusters. Redshift evolution of these trends will be 
investigated in future work. 

This paper is organised as follows. In Section [5] we 
briefly describe the N-body simulations we use and how we 
compute a number of halo properties. The method for the 
identification of the halo environment is presented in Sec- 
tion [3] together with a number of tests that show how well 
the method performs. Our main results on the environmen- 
tal dependence of the halo properties are given in Section U 
Finally, Section [5] summarises the main conclusions that we 
draw from our work. 



3.8 x 10 8 /i _1 M o and 3.0 x 10 9 /i _1 A'/ Q , respectively. The 
simulations follow the evolution of Gaussian density fluctu- 
ations characterised by a scale-free initial power spectrum 
with spectral index n — 1 and normalisation as = 0.9 (with 
as the rms linear density fluctuation within a sphere of 8 ft" 1 
Mpc comoving radius). The ini tial conditions wer e gener- 
ated using the GRAFIC2 tool (|Bertschingerll200ll ) for the 
redshift z at which the rms density fluctuation on the small- 
est resolvable scale in each box equals 0.1. This corresponds 
to z ~ 79, 65 and 52 for L\ , L2 and L3 respectively. Par- 
ticle positions and velocities were saved for 30 time-steps 
logarithmically spaced in expansion parameter a between 
z = 10 and 2 = 0. 



2.1 Halo identification and properties 

Virialised dark-matter haloes were identified using the stan- 
dard friends-of-friends (FOF) algorithm with a linking 
length equal to 0.2 times the mean inter-particle distance. 
We only considered haloes containing at least 300 particles, 
since virtually all of the halo properties we investigated show 
strong numerical artefacts when measured for less well re- 
solved haloes. We found 13353, 16296 and 21041 of such 
haloes for L\, L2 and L3, respectively. The most massive 
groups in the three simulations contain nearly 10 6-7 parti- 
cles and have masses 4.3 x lQ 1Jk h~ x M & , 7.6 x 10 14 /i _1 M© 
and 2.2 x 10 15 /i _1 Mq. Our catalog therefore spans five or- 
ders of magnitude in halo mass with high resolution haloes, 
ranging from the size of dwarf galaxies to massive clusters. 

We characterised the mass assembly and merging his- 
tory of the halos as follows. For each halo at redshift 2, we 
identified a progenitor at z p > 2 by intersecting the sets of 
their particles. The main progenitor was then chosen to be 
the most massive halo at each redshift that contributes at 
least 50 per cent of its particles to the final halo. We then 
defined the formation redshift 2f orm as the epoch at which a 
main progenitor which has at least half of the final mass first 
appears in the simulation and interpolated linearly between 
simulation snapshots in log 2 to find the point where exactly 
half of the mass is accumulated. 



2.1.1 Halo Shapes 

In order to quantify the shape of FOF haloes, we determined 
their moment of inertia tensor, defined as 



Ijk = m 



jk %i,j%i,k ) 



(i) 



2 N-BODY SIMULATIONS 



We used the tree-PM code GADGET-2 <|Springelll2005r i to 
follow the formation and the evolution of the large-scale 
structure in a flat ACDM cosmology. We have assumed the 
matter density parameter Q m = 0.25, with a baryonic contri- 
bution fib = 0.045, and the present-day value for the Hubble 
constant Ho = 100/t km s _1 Mpc -1 , with h = 0.73. In par- 
ticular, we performed three N-body simulations, each con- 
taining 512 3 dark matter particles in periodic boxes of size 
Li = 45 ft" 1 Mpc, L 2 = 90 h' 1 Mpc and L 3 = 180 ft" 1 Mpc. 
The corresponding particle masses are 4.7 x 10 7 /i _1 M Q , 



where m is the particle mass, n = (xi,i, Xi t %, Xi,z) is the dis- 
tance of the i-th particle from the centre of mass of the halo 
and 5jk denotes the Kronecker symbol. The eigenvectors of 
I are related to t he lengths of the principal axes of inertia 
h > h > h (e.g. IBett et~aill2006h . We used the following 
dimensionless quantities 



h 



and 



T = 



il_-H 



'3 



(2) 



to measure sphericity and triaxialit y of the haloes (e.g. 
iFranx et al"1ll99ll ; IWarren et al.lll992T ). A spherical halo has 
S = 1, a needle S = 0, a prolate halo T = 1 and an oblate 
one T = 0. 
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2.1.2 Halo Spin Parameter 

The spin par ameter o f a ha lo is a dimensionless quantity 
introduced bv lPeeblesI (|l969l ) that indicates the amount of 
ordered rotation compared to the internal random motions. 
For a halo of mass M and angular momentum J it is defined 

as 



1/2 



A 



I J I \E\ 
CM 5 /? 



(3) 



where the total energy E — T + U with T the kinetic en- 
ergy of the halo after subtracting its bulk motion and U the 
potential energy of the halo produced by its own mass distri- 
bution. Determining the potenti al energy of massive haloes 
is computationally expensive, so IBullock et all (|200lf) intro- 
duced the alternative spin parameter 



A' 



I J vir 



~\/2iM v i T Vvir-Rvir 



(4) 



Here all quantities with the subscript "vir" (angular mo- 
mentum, mass and circular velocity) are computed within a 
sphere of radius i? v ir which approximates the virial radius 
of the halo. As this quantity is not well defined for FOF 
groups, we took R v i r to be a fraction a of the maximum 
distance between a halo particle and the centre of mass. To 
accommodate possible fuzzy boundaries of the haloes, we 
chose a value of a = 0.95. We verified that the particular 
choice of a does not have an impact on the distribution of A' 
and rem ains unchanged even when going as low as a = 0.1 
(see also IBullock et all [ioOlh . Under the assumption that 
the halo is in dynamical equilibrium, V^ r = GM V1T / R v i r , 
the spin parameter can be rewritten as 



A' 



J vir 



3/2 ' 



(5) 



We found a spurious increase in A' for haloes consisting of 
less than 250-300 particles. This numerical effect occurred 
for all of our three simulated boxes. The median spin A^ed is 
roughly 10 per cent higher for haloes with only 100 particles 
than for haloes consisting of more than 300 particles. 



3 ORBIT STABILITY AND ENVIRONMENT 
3.1 Basic theory 

We use a simple stability criterion from the theory of dy- 
namical systems to distinguish between haloes residing in 
clusters, filaments, sheets or voids. Consider a test particle 
moving in the peculiar gravitational potential, <j>, generated 
by a cosmological matter distribution frozen in time (e.g. no 
Hubble drag). The equation of motion in comoving coordi- 
nates for this test particle is x = — V(j>, where the dot repre- 
sents derivatives with respect to a fictitious time. Assuming 
that at the centre of mass of each halo x, the gravitational 
potential has a local extremum (i.e. V0(xj) = 0), the fixed 
points of the test particle equation of motion are exactly at 
the points x^. We can thus linearise the equation of motion 
at the points ij and find the linear system 



where the tidal field Ty is given by the Hessian of the grav- 
itational potential 



didj 4> ■ 



(7) 



Thus the linear dynamics near local extrema of the gravita- 
tional potential is fully governed by the three (purely real, as 
Tij is symmetric) eigenvalues of the tidal field tensor. We use 
the number of positive eigenvalues of Tj to classify the four 
possible environments a halo may reside in. Note that the 
number of positive eigenvalues is equivalent to the dimen- 
sion of the stable mani fold at the fixed points. In analogy 
with Zel'dovich theory (|Zel'dovichlll97Ch . we define as 

(i) voids the region of space where Tij has no positive 
eigenvalues (unstable orbits); 

(ii) sheets the set of points with one positive and two neg- 
ative eigenvalues (1-dimensional stable manifold); 

(iii) filaments the sites with two positive and one negative 
eigenvalue (2-dimensional stable manifold); 

(iv) clusters the zones with three positive eigenvalues (at- 
tractive fixed points). 

Dropping the assumption of local extrema of the gravita- 
tional potential at the centres of mass of the haloes intro- 
duces a constant acceleration term to the linearised equa- 
tions of motion. This zeroth-order effect can be dispersed of 
by changing to free-falling coordinates. The deformation be- 
haviour introduced by the first-order term, however, remains 
unchanged. 



3.2 Implementation 

In order to determine the eigenspace structure of the tidal 
field tensor, we first compute the peculiar gravitational po- 
tential <f> from the matter density distribution via Poisson's 
equation 



= 4-KGpS, 



(8) 



where p and 5 respectively denote the mean mass density 
of the universe and the overdensity field. For our N-body 
simulations, we solve Poisson's equation using a fast Fourier 
transform on a grid of twice the particle resolution (1024 3 
grid cells). The density field S is obtained by using Cloud- 
In-Cell interpolation of the particles onto the grid and then 
smoothed using a Gaussian kernel Kr s . In this case, the 
smoothing length, R s , and the mean mass contained in the 
filter, M s , follow the relation 



Rs 



P 



1/3 



(9) 



To solve Poisson's equation on the grid, we apply the Green's 
function of the symmetric 5-point finite difference op- 
erator that we later use to compute the tidal tensor. Alto- 
gether, we hence find the solution for the smoothed gravita- 
tional potential through the double convolution 



S * Kr b * G 



(2) 



(10) 



-Ty(Xfe) {Xj - X k ,j) ■ 



(6) 



We then apply the second derivative operator to (j>R s and get 
the diagonal components of the tidal tensor. For the off-trace 
components, we apply twice the symmetric first derivative 
operator in the corresponding coordinates. Although the sec- 
ond derivative operator cannot be produced from applying 
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Figure 1. The volume fraction being classified as clusters, fila- 
ments, sheets or voids for our 180 h' 1 Mpc box as a function of 
the smoothing scale R B . The vertical dotted line at R B = 2.1h~ 
Mpc indicates the smoothing scale adopted in this paper. The 
solid grey lines indic ate the predicted vo lume fractions for a Gaus- 
sian random field jPoroshkevichl Il970h . For very large J2 S , the 
non-Gaussian density field of the simulations asymptotes to the 
predicted fractions of 42 per cent for sheets and filaments and 8 
per cent for voids and clusters. Volume fractions are evaluated on 
a 128 3 Cartesian subgrid. 

twice a symmetric first derivative operator, following this 
scheme ensures that the trace coincides with the smoothed 
overdensity to machine accuracy, while the off-trace com- 
ponents are indeed symmetric and are not suffering from a 
spurious self-potential. Finally, we compute the eigenspace 
structure of the tensor at each halo's centre of mass. 

3.3 Optimisation 

Our criterion for determining the halo environment contains 
one free parameter, namely the smoothing radius of the 
Gaussian kernel, R s . This corresponds to the typical length- 
scale over which we determine the dynamical stability of the 
orbits. The particular choice of R s directly affects the local 
eigenstructure and thus changes the classification of environ- 
ment. Smoothing on the scale of single haloes picks out each 
single halo as a stable cluster in the sense of the definition. In 
Figure Q] we show how the choice of R s affects the fraction of 
the simulated volume classified in the four categories. These 
fractions continuously vary with R s , which implies that some 
haloes change their classification. For R B 3> 10/i -1 Mpc, the 
density field becomes approximately Gaussian, and we ob- 
serve convergence to the theoretical volume fractions of 42 
per cent for sheets and filaments, and 8 per cent for voids 
and clusters (Doroshkevich 1970, see also Shen et al. 2006). 
To illustrate the transition of haloes between environment 
classes, in Figure [2] we use a snapshot at z — to high- 
light the haloes that are assigned to different environments 



when the smoothing scale is changed from 2.1 h~ x Mpc to 
4.5 hT 1 Mpc (corresponding to a change by a factor of 10 in 
M s ). Basically, increasing the smoothing scale: i) increases 
the number of haloes in voids at the expenses of the sur- 
rounding sheets (panel b); ii) moves the thin filaments sur- 
rounding a thicker one from the sheet environment to the 
filament one (panel c); in) moves the thin filaments sur- 
rounding a void from the filament environment to the sheet 
one (panel d); iv) increases the size of massive clusters lo- 
cated at the intersection of filaments at the expenses of the 
ending points of filaments themselves (panel e); v) moves 
the densest clumps located along filaments from the clus- 
ter environment to the filament one (panel f). Table [T] lists 
the fraction of the total number of haloes that are assigned 
to the 16 possible classifications with the two smoothing 
scales. The halos that contribute to the off-diagonal ele- 
ments of this "transition matrix" typically live in regions 
where the tidal field has one nearly vanishing eigenvalue. 
In these transition regions, a modification of Rs can easily 
change the sign of this eigenvalue and thus the association 
of the corresponding halo to its environment. This results 
from using sharp boundaries (positive vs negative eigenval- 
ues) to classify the different environments. Note that only 
a negligible fraction of the haloes inverts the sign of more 
than one eigenvalue of the tidal field when the smoothing 
scale is changed, indicating that our classification is indeed 
physical. Based on Figure [2] we conclude that the combined 
use of two (or more) smoothing scales can be used to clas- 
sify a larger variety of environments with respect to the 
basic four that can be found with a fixed resolution, and 
in particular to identify boundary regions that bridge be- 
tween the basic four types. We will explore this potentiality 
of the orbit-stability method in future works. For simplic- 
ity, in this paper we only consider a single smoothing scale, 
Rs = 2.1 h" 1 Mpc (corresponding to M B « 1Q 13 /i _1 M ) 
which provides startling agreement between the outcome of 
the orbit-stability criterion and a visual classification of the 
large-scale structure. The resulting classification of halo en- 
vironments is highlighted in the top-left panel of Figure [2] 
using different colours. For R B = 2.1 Mpc, the volume 
fractions occupied by voids, sheets, filaments and clusters 
are, respectively, 13.5%, 53.6%, 31.2% and 1.7%. This sug- 
gests that we identify as voids just the inner parts of the 
most under-dense regions (see also Figure [3} and consider 
as sheets the volume-filling regions around them. At the 
same time, our clusters always contain haloes with a virial 
mass M v i r > 10 13 /i _1 Mq and, in some cases, haloes with 
M vir > lO 14 /i _1 Af and radius R vil > lh' 1 Mpc (which 
are usually tagged as clusters). These haloes typically con- 
stitute the central parts of what we identify as clusters. By 
definition, our "cluster environment" extends to distances 
which are significantly larger than _R v ir and also includes all 
the smaller haloes that are infalling onto or orbiting around 
the central one. For the value of R s adopted in this paper, 
we find that "our" clusters have typical diameters of a few 
Mpc. We have tested that all our findings do not depend on 
the precise choice of R s . 

3.4 Orbit stability vs density 

Most of the work on the environmental dependence of halo 
properties has hitherto considered the local density as a 
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Figure 2. Classification of halo environments in a slice of 10/i -1 Mpc thickness for the 180/i _1 Mpc box. Panel a) shows the classification 
scheme for a smoothing scale of R B = 2.1h~ x Mpc with the following colour coding: clusters (red), filaments (blue), sheets (green) and 
voids (orange). Panels b) to f) show in red those haloes that change classification in a specific way when the smoothing scale is increased 
to R s = 4.5/i -1 Mpc, all other haloes are represented in grey. Panel b) represents sheets at smaller smoothing R B = 2.1/i — 1 Mpc that 
become voids at larger smoothing R B = 4.5?t — 1 Mpc. Panel c) shows sheets that become filaments, panel d) filaments that become sheets, 
panel e) filaments that become clusters and panel f) shows clusters that become filaments. To achieve higher spatial resolution in the 
visual representations, all haloes down to 10 particles are shown in the panels above. 
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Table 1. Transition matrix for halo classification between 
smoothing at M a = 1O 13 M0, indicated by "(S)", and M s = 
lO 14 Af0, indicated by "(L)". Matrix entries are given in per cent 
of the total number of haloes. Non-diagonal elements represent 
haloes that change classification. 
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Figure 4. Mass function of the haloes residing in voids, sheets, 
filaments and clusters. Abundances in the whole box have been 
rescaled by the corresponding volume fractions. The solid grey 
line represents the total mass function, not split into environ- 
ments. Haloes from all three simulations are included. The total 
mass function perfe ctly coincides with common analytic fits (e.g. 
I Jenkins et al]|200l)V 



4 HALO PROPERTIES AND ENVIRONMENT 

In this section we present a detailed study of halo properties 
at 2 = as a function of the cluster, filament, sheet and void 
environment determined by our orbit stability criterion. 



Figure 3. Volume weighted probability distribution of the local 
density for clusters, filaments, sheets and voids. Statistics have 
been obtained combining all three simulation volumes. Statis- 
tics weighted by halo abundance shifts the distributions to over- 
densities roughly a factor of 2 higher. Note that the stability 
criterion naturally finds "clusters" in the highest density regions 
and "voids" in the lowest and thus disambiguates any definition 
of environment that is solely based on density measures. 



measure of environment (e.g. Lemson fc Kauffmannl 1 19991 : 
IMaccio etal]|2006l : iMaulbetsch et al.ll2006h . Density corre- 
sponds to the trace of the tidal field tensor Tij, and thus 
provides more limited information regarding the dynamical 
properties of the local flow compared to our classification, 
which is based on all three eigenvalues. In Fig. [3] we show 
that local overdensity is largely degenerate relative to the 
four categories we derive from the eigenstructure. Density 
correlates with the dimension of the stable manifold, e.g. 
the median overdensity in each environment is -0.79, -0.55, 
0.28 and 4.44 for voids, sheets, filaments and clusters, re- 
spectively. However, it is not possible to recover, from the 
density field, the more detailed environmental information 
that we derive from the tidal field tensor. A simple envi- 
ronmental classification that is based on density therefore 
mixes our halo populations. 



4.1 Mass function 

Figure [4] presents the mass functions of the haloes resid- 
ing in the different environments. The low-mass end has the 
same slope in all environments, but the position of the high- 
mass cutoff is a strong function of environment. The cluster 
mass function is top-heavy with respect to voids, while fil- 
aments and sheets lie in between. As expected, the mean 
halo density is higher in clusters and lower in voids. All this 
is in good qualitative agreement with the conditional mass 
function as a function of local density derived from analytic 
models l|Bond et al.lll99ll : lBower]|l99ll ). 



4.2 Halo Shapes 

Figure [5] shows the median of the shape parameters S and 
T for haloes in filaments and clusters as a function of their 
mass (the void sample contains too few haloes and the 
sample for sheets shows an identical behaviour to the fil- 
aments). Only when requiring that haloes in the samples 
contain at least 500 particles we find convergence of the me- 
dian shape parameters at the lower mass end. The overall 
mass dependence of S and T is in good agreement with 



iBett et all 120061 : IMaccio et al.1 12006ft . lAllgood et all (|2006l ) 



previous st udies (e.g. lAllgood et al.ll2q06j : lAltav et a l1l2006l : 



fit a power law to the median S as a function of halo mass, 
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Figure 5. Median halo sphericity (left) and triaxiality (right) as a function of halo mass for haloes in filaments and clusters. The 
behaviour for haloes in sheets is almost identical to that for filaments. The shaded area ind icates the central la scatter in the whole 
sample, not split by environment. The dark grey lines indicate the fits of lAlleood et al" (2006) for S and lBett et al. (2006) for S and T, 
the black lines show our fits to haloes with masses M < 2 X W 12 h~ 1 Mq. Parameters are given in section [4.21 
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Figure 6. Median mass of the main progenitor of haloes in the 
mass range 5 X lO lo /i -1 M < M < 5 X 10 11 h-lM Q over red- 
shift as a fraction of the mass at z = for the four different en- 
vironments. The shaded area indicates the la spread of haloes in 
filament environments. The spread is slightly larger for haloes in 
clusters. The dotted grey line indicates where Zf ornl is measured. 



well with the mea sured slope s of bo th fitting formulas for 
masses M > M c . iBett et all (l2006h argue th at the offset 
of their fit with respect to Allgood et al.l l|2006t ) results from 
different halo finding algorithms which also explains why our 
haloes are slightly less spherical. We do not find evidence for 
decreasing spher icity at the low-mass end as indicated by 
I Bett et al . ( 2006), based on haloes with less than 300 parti- 
cles. However, we clearly detect a decrease in slope for halo 
m asses M < 10 12 /i~ 1 M b with respect to the fitting formula 
of I Allgood etafl (I2006T I. 

The vast dynamic range of our suite of simulations allows the 
unprecedented exploration of the low-mass end with high- 
resolution haloes (500 particles for a 2 x 10 10 Ti -1 Mq halo). 
For masses in the range 2 x lO 10 /! -1 !^ < M < M c we de- 
tect a clear dependence on environment. Haloes residing in 
clusters tend to be less spherical and more prolate almost 
independently of mass. In contrast, haloes in filaments tend 
to be slightly more oblate as one might expect from accre- 
tion of matter onto the filament. The difference between the 
two classes are, however, small with respect to the intrinsic 
scatter. For masses M < M c , our results are well described 
by a fit of the following form: 



Smcd 



S2 l 
Si + lo; 



M 



= h + 



100 610 lO^/i-iMg 
fa , M 
100 ° Sl ° lO^ft-iMn ' 



(11) 
(12) 



while IBett et al.1 (|2006h detect a change in slope at masses 

M c w 2 x 10 12 /i _1 M Q . This breakpoint M c is present also in Sl = 0.66 ± 0.08, 

our findings. Interestingly, it coincides with the mass above s 2 = —3.6 ± 0.7; 

which we do not find any significant dependence of the ti — 0.66 ± 0.03, 

shape parameters on environment. Our results agree very ti = 2.4 ±0.24 



in filaments, M < M c , (13) 
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Figure 7. The median formation redshift z mc( j for haloes from 
our three simulated boxes as a function of their mass. Errorbars 
indicate the error in the median. The grey line indicates the re- 
sult of a robust fit to the displayed medians. For masses below 
5 X 10 12 h~ 1 Mq we find a strong correlation with our definition of 
environment. The black lines indicate robust fits to the values for 
haloes with masses M < 5 X 10 12 ?i — 1 Mq that reside in the corre- 
sponding environments. The fit parameters for all environments 
are given in section [4.31 



Si 
S2 
ti 

t 2 



0.64 ± 0.05, 
-1.5 ±0.42; 
0.69 ±0.03, 
1.0 ±0.24 



in clusters, M < M c , (14) 



where M c = 2 x 10 12 /i _1 Mq. These values are obtained with 
a robust iterative least-squares fit using a bisquare estima- 
tor. 



4.3 Assembly history and formation redshift 

In Figure [6] we show the assembly history of haloes with 
masses 5 x lO lo /i _1 M < M < 5 x 1O 11 /i" 1 M in the differ- 
ent environments. In particular, we plot the median mass of 
the main progenitor as a function of redshift at which it is 
identified. The shaded area indicates the central la spread 
for haloes in the filament environment. Although haloes tend 
to assemble their mass earlier in clusters and later in voids, 
the effect is relatively small with respect to the intrinsic 
sc atter. This is in very go od agreement with the findings 
of iMaulbetsch et al.l (|2006l ). These authors investigated the 
mass assembly history splitting the halo sample by density, 
smoothed on 4/i _1 Mpc. Their high density sample (S > 5) 
roughly corresponds to our densest clusters, while the low 
density sample (5 < 0) includes voids, sheets and the lower 
density filaments. 

Both shape and scatter of the mass assembly curve de- 
pend strongly on the mass range at which they are evalu- 
ated. This points to a strong relation between the formation 
redshift of haloes, their mass and environment. In Figure [7] 



Figure 8. The distribution of formation redshifts for haloes with 
masses 2 X 10 10 h- 1 M Q < M < 10 11 /i _1 M Q split into our four 
environment categories. Note that there are many more haloes in 
filaments than in clusters in this mass range. 



we plot the median formation redshift z mc d as a function of 
halo mass M for the haloes from our three simulations. Er- 
rorbars are estimates of the error in the median computed 

as 



ZQ.84 — Z0.16 



(15) 



where zo.84 and zo.w denote the 84th and 16th percentile 
of the distribution of Zf orm , corresponding to the la spread 
if the underlying distribution were Gaussian, and Ah is the 
number of haloes used to sample the distribution. For four 
decades in mass, ranging from 10 10 /i -1 Mq to 10 14 /i -1 Mq, 
we find a tight logarithmic relation between halo mass and 
formation redshift, reflecting the hierarchical structure for- 
mation paradigm. Results for all boxes agree very well when 
considering haloes of at least 300 particles. We fit a function 
of the form 



Zmcd — Cl 



C2 log'i 



M 



lO 12 /;,- 1 Mi 







(16) 



The parameters given by a robust fit to all haloes from the 
three simulations are: 

ci = 1.29 ±0.07, 
c 2 = 0.312 ± 0.006. 

For haloes with masses between 10 10 /i _1 Mq and w 
10 12 /i _1 Mq we find that z mc d strongly depends on envi- 
ronment. This dependence increases, the lower the mass of 
the haloes. Our r esults a re in very g ood a g reement with 
ISheth fc Tormenl (1200-d ) iGao et all l|2005h . lHarker et aD 
l|2006h and iReed et all (120061 '). These authors found that 
haloes of given mass but different formation epochs show dif- 
ferent clustering properties. In particular, they have shown 
that small-mass haloes with higher formation times cluster 
more strongly and are thus most likely associated to denser 
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environments. For haloes with masses M < 5 x 10 12 /i _1 Mq 
we again fitted relation (|16|) separately for cluster, filament, 
sheet and void environments. The slope parameters C2 are 
significantly different for the four environments. A robust fit 
to the data combined from all three simulations yields the 
fit parameters 



clusters, M < 5 x 1O 12 /i _1 M ; 



Cl 


= 1.42 ±0.39 


C2 


= 0.54 ±0.03 


Cl 


= 1.30 ±0.04 


C2 


= 0.39 ±0.01 


Cl 


= 1.21 ±0.11 


C2 


= 0.28 ±0.01 


Cl 


= 1.36 ±0.48 


C2 


= 0.08 ±0.04 



filaments, M < 5 x 1O 12 /i _1 M ; 
sheets, M < 5 x W 12 h' 1 M Q - 
voids, M < 5 x 10 u ft _1 Mo. 



For M > 5xlO 12 /i _1 M , we do not find any dependence on 
environment and the relation between z me d and halo mass 
is best fit by the relation for all haloes given above. 

The differences between the environments become even 
more significant when considering mean values of Zf orm in- 
stead of the medians due to the skewness of the formation 
redshift distributions in each mass bin. To illustrate this, we 
plotted in figure [8] the distribution of formation redshifts for 
haloes with masses 2 x W 10 h~ 1 M Q < M < 10 u /i -1 Af Q 
in the four enviro nments. In very good agreement with 
IWang et all (|2006l ). we find that the oldest haloes with 
£form > 3 are relatively overrepresented in cluster envi- 
ronments. Thus, small mass haloes in the vicinity of clus- 
ters tend to be older, and there has to be some effect that 
prevents them from strong continuous accretion and major 
mergers in this environment. In absolute numbers, we find 
a comparable amount of these very old low mass haloes also 
in our filament environments, such that, to a l esser extent, 
a sim ilar effect must be present in filaments. IWang et al.l 
(2006) suggest that the survival of these fossil haloes may 
be related to the "temperature" of the surrounding flow. 
It is evident from these findings, that this "temperature" 
would then strongly correlate with the dimension of the sta- 
ble manifold in our classification of environment. The pre- 
cise connection has to be investigated in future work, but it 
is conceivable that the higher the number of stable dimen- 
sions, the less coherent and more accelerated is the infall of 
surrounding matter, and the stronger the heating of dark 
matter random motion. 



4.4 Halo Spin 

We investigate the dependence of the halo spin parameter 
A' on environment. Figure [9] shows the distribution of A' in 
the mass ranges 5 x 10 10 h~ 1 M Q < M < 5 x 10 11 /i" 1 M Q and 
M > 5 x 10 12 ft _1 M©. In the high -mass bin, the distribution 
of spin parameters is well approximated by a log-normal 
probability density function, 



p(A') = 



1 



A' ay V27T 



exp 



log 2 (A7A&) 
2 a 2 



(17) 



with best-fitting parameters Aq = 0.035 and width ay = 
0.70. However, for M < 5 x 10 12 /i _1 M Q , we find a tail of 
rapidly spinning haloes that is most prominent in clusters, 
and to a lesser extent in filaments. For the mass range 5 x 
W 10 h~ 1 M Q < M < 5x W 11 h~ 1 M Q , we find good agreement 



of all environments with a log-normal distribution only for 
spin parameters A' < 0.1. The fit parameters for A' < 0.1 
in the low mass regime are X' — 0.030 and ay = 0.61. 
Our findings for the parameter An agree well wi th earlier 
findings (e.g. iBullock et alj|200ll ; iBett et al.ll2006l ). At A' « 
0.1, however, we detect evidence for a departure from the 
log-normal distribution that is very well fit by a power-law 
behaviour. We find 



p(X'\X' > 0.1) = 0.0012 A' ~ 31 
for haloes in filaments and 
p(A'|A' > 0.1) = 0.035 A'" 18 



(18) 



(19) 



for haloes in clusters. This tail is almost independent of 
the assumed value of a, i.e. the fraction of the virial radius 
within which A' is determined. However, the environmental 
dependence of the spin distribution slightly decreases when 
only the very innermost parts of a halo are used to determine 
A'. We have also verified that the high-spin tail of the dis- 
tribution is not affected by measurement errors of the halo 
spin, i.e. the statistics remains unaltered when only haloes 
containing > 1000 particles are considered. 

Our results appea r to be in disagreement with 
lAvila-Reese et all (2005) who found that haloes in clus- 
ters are less rapidly spinning than in the field. However, 
a direct comparison is problematic since a) we use a differ- 
ent halo finder algorithm, b) we do not consider sub-haloes 
(which likely suffer strong tidal stripping), and c) we use 
a different definition of the cluster environment. However, 
we agree well with their finding that the parameter ay of 
the log-normal fit is significantly larger for haloes in clus- 
ter environments than for h aloes in under-dense regions. 
iHetznecker fc BurkertJ (2006) have recently shown that the 
halo spin parameter increases significantly after a major 
merger, and relaxes to more standard values after 1-2 Gyr. 
The formation redshift of a halo, as defined in section [5] 
is a good indicator for the occurrence of major mergers. 
Low formation redshifts correspond to recent major merg- 
ers, while high values of Zf om denote less violent accretion 
histories. Figure [10] shows the median spin parameter as a 
function of z Iorm for two mass-bins. We find that, in all en- 
vironments and mass ranges, the median A' is a decreasing 
function of z{ olnl . At the same time, for a given Zform, the spin 
parameter shows an important environmental dependence: 
low-mass haloes tend to spin faster if they reside in clusters 
while massive haloes tend to spin slower in this environ- 
ment. The haloes with the largest spin parameter (median 
A' > 0.1) are low mass haloes M < 5 x 1O 12 /i _1 M that re- 
side in clusters and have zf olm < 1. However, for fixed Zf orm 
haloes in clusters have higher median A' compared to the 
other environments. 



4.5 Angular Momentum Alignments 

Do halo spin directions retain memory of the cosmic web in 
which the haloes formed? Both filaments and sheets have a 
preferred direction given by the structure of the eigenspace. 
While filaments are one-dimensional structures with a pre- 
ferred direction in space, sheets are two-dimensional and can 
thus be uniquely described by their normal vectors. Using 
the definition in section 13.11 these directions are given by 
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Figure 9. Left panel: Distribution of halo spin parameter A' for haloes in the mass interval 5 x 10 10 /i _1 Af Q < M < 5 x I0 11 h- 1 M Q 
residing in clusters, filaments, sheets and voids. Statistics are combined for all three simulation volumes. The solid grey line indicates 
the fit of a log-normal distribution to the sample A' < 0.1, not split by environment. The dashed grey line shows a power-law fit to 
the distribution of A' > 0.1 for haloes in filaments, the dash-dotted grey line the corresponding fit to cluster haloes. Right panel: Spin 
parameter distribution for halo masses M > 5 X 10 12 ?t — 1 Mq for which haloes in voids are not present. The solid grey line shows the fit 
of a log-normal distribution to the whole sample, not divided into environment types. All fit parameters are given in section [4.41 
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Figure 10. The median spin parameter A' of haloes in the mass range 5 X 10 10 /i _1 Mq < M < 5 X lO 11 ^ - 1 Mq (left) and M > 
5 X 10 12 h~ 1 Mq (right) for the four different environments as a function of their formation redshift. Errorbars indicate the lc uncertainty 
in the median. 



the unit eigenvector v corresponding to the negative eigen- 
value of the tidal field tensor for filaments and the positive 
eigenvalue for sheets. One can therefore compute the de- 
gree of alignment between the angular momentum vector 
of a halo and the respective eigenvectors of the environ- 
ment in which it resides, J • v. Figure [11] shows the dis- 



tribution of alignments between halo angular momentum 
and both filament direction and sheet normal vector in the 
two mass-bins 5 x 10 10 /i _1 Mq < M < 10 12 h~ 1 M (s and 
M > 1O 12 /i _1 M0. For haloes in filaments we find only a 
weak trend for their angular momenta to be aligned with the 
filament direction. Haloes in sheets, however, show a very 
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Figure 11. The alignment between halo angular momentum vec- 
tors and the eigenvector corresponding a direction perpendicu- 
lar to the sheets (top), and corresponding the direction of the 
filaments (bottom), for haloes in these two environments. Halo 
populations are divided in two bins 5 X 10 10 h~ 1 Mq < M < 
10 12 h- 1 M e and M > 10 12 h- 1 M Q . The dotted grey lines indi- 
cate a random signal. 



strong tendency to have their angular momentum parallel 
to the sheet. Si milar correlations are also found in wa lls de- 
limiting voids (jPatiri et all 120061 : iBrunino et alj 120061 ). and 
might be reflected i n the distribution of galactic disks (e.g. 
iNavarro et all 12004 iTruiillo et all |2006j ). We did not de- 
tect any strong correlation with eigenvectors of the other 
environments. The presence of alignments between large- 
scale structures and halo spins could produce a coherent 
alignment of galaxy shapes and thus generate a systematic 
contamination in weak lensing maps of co smic shear (e.g. 
iHirata fe SeliaJj|2004 iHevmans etal]|2006h . 

We next compute correlations of the intrinsic angular 
momentum of each halo with both the intrinsic angular 
momentum and the orbital angular momentum of neigh- 
bouring haloes residing in the same e nvironment. We define 
the spin-spin correlation function as (|Porciani et al . 2002a; 
iBailin fc Steinmet dliooBh : 



£j.j(r) = (J(x).J(x + r)>, 



(20) 



where J is the intrinsic angular momentum of each halo, 
and the average is taken over all pairs of haloes which are 
separated by a distance r and reside in the same environment 
class. Similarly, we define the spin-orbit correlation as 



£j. L (r) = <J(x)-L(x + r)>, 



(21) 



where L is the relative orbital angular momentum be- 
tween two haloes separated by a distance r. Figure 1121 
shows the spin-spin correlation for haloes in two mass bins, 
5 x lO lo ft _1 M < M < 5 x IO^/i^A/q and M > 5 x 
1O 12 /i _1 M . We find a significant correlation only for haloes 
with M > 5 x 10 12 h~ x Mq in cluster environments. These 
haloes have a strong tendency to have their spin vectors 



Figure 12. The mean alignment of intrinsic spin angular mo- 
menta between haloes in filaments, panels a) and c), and clus- 
ters, panels b) and d). Data is plotted for the two mass- bins 
5 X 10 10 h- 1 M e < M < 5 X W^h^M®, panels a) and b), and 
M > 5 X 10 12 h~ 1 Mq, panels c) and d). Errorbars are the la 
uncertainty of the mean. The grey line indicates the mean corre- 
lation for the whole halo population, independent of environment, 
for the same mass bins. The dotted line represents a random sig- 
nal with no correlation. 



anti-parallel to the spins of haloes within a distance of a 
few Mpc. All other correlations are essentially consistent 
with a random signal. The results for haloes with masses 
5 x IO^^Mq < M < 5 x 10 12 ft~ 1 M Q are fully consistent 
with those for the lower mass bin and therefore not shown 
in the plots. Regarding the alignment of spin and orbital an- 
gular momenta, the results, given in Figure [131 show a much 
stronger signal and a clear dependence on environment. We 
find an evident tendency for the two angular momenta to 
be parallel regardless of mass and environment. Remarkably 
this correlation significantly extends out to ~ 2/i _1 Mpc in 
all environments and is most prominent for smaller haloes 
in clusters and massive haloes in filaments. 



5 SUMMARY 

We have presented a new method to classify dark-matter 
haloes as belonging to four different environments: clusters, 
filaments, sheets and voids. This scheme computes the di- 
mensionality of the stable manifold for the orbits of test 
particles by simply looking at the number of positive eigen- 
values of the local tidal tensor. The algorithm contains only 
one free parameter: the smoothing radius for the gravita- 
tional potential. This quantity fixes the length-scale over 
which the stability of structures is determined and can be 
fine tuned to optimise the classification. At the same time, 
combining the results obtained adopting two or more differ- 
ent smoothing scales allows us to select regions with partic- 
ular properties in the large-scale structure (e.g. transition 
regions between the basic four environments). 
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Figure 13. The mean alignment of intrinsic spin and relative 
orbital angular momentum between haloes in filaments, panels a) 
and c), and clusters, panels b) and d). Data is plotted for the two 
mass-bins 5 X 10 w h~ 1 M e < M < 5 X 10 n h~ 1 M @ , panels a) and 
b), and M > 5 X lO 12 h. _1 M0, panels c) and d). Errorbars are 
the lex uncertainties of the mean. The solid grey line indicates the 
mean correlation for the whole halo population, independent of 
environment, for the same mass bins. The dotted line represents 
a random signal with no correlation. 

Our classification scheme correlates with local density 
so that the densest regions are always associated with clus- 
ters and the emptiest with voids. However, our method re- 
tains more information on the local dynamics and a simple 
halo classification based on density will unavoidably mix our 
populations up. 

We have used the classification scheme to study how 
the properties of isolated dark-matter haloes depend on the 
environment in which they reside at z — 0. Our main results 
can be summarised as follows. 

1) Halo shapes 

• Massive haloes with M > a few x 10 12 Mq do not 
show any significant dependence of their shape on environ- 
ment. 

• Less massive haloes in clusters are less spherical and 
more oblate than in other regions but the trend is generally 
weak compared with the intrinsic scatter. 

2) Halo formation times 

• For the whole halo population (not split by environ- 
ment) we found a very strong correlation between median 
formation redshift and halo mass. A fit to this relation 
which holds for halo masses between lO lo /i _1 M < M < 

is given in equation (|16p . This dependence is a 
direct consequence of hierarchical structure formation. 

• For M < 5 x 10 12 /i _1 M Q haloes of fixed mass in the 
four environments have significantly different mass assem- 
bly histories. In particular, cluster haloes tend to be older 
while void haloes younger. All this hints at mechanisms that 



suppress the growth of lower mass haloes in clusters and 
lead to an enhan ced survival rate of fossil haloes (see e.g. 
Wan g et al . 2006, for a possible explanation). 

• Analytic fitting formulae for the dependence of the me- 
dian formation redshift on halo mass and environment are 
given in Section [4.31 

3) Halo spins 

• The median spin parameter of all haloes is the highest in 
clusters followed in order by filaments, sheets and voids. This 
dependence, presumably, has its origin in the tidal-torque 
history of the haloes which likely correlates with the specific 
eigenstructure of the tidal field at the final halo position 
ijBond et alJll99&|Porciani et al.ll2002al ibh. 

• This trend is reversed for massive objects. Haloes with 
M > 5 x 10 12 /i _1 Mq in clusters are less rapidly spinning 
than in filaments. 

• On the other hand, for smaller masses, haloes in clusters 
generally possess higher spin parameters than in the other 
three environments. As these rapidly spinning haloes have 
also the most recent formation time, we conjecture that the 
high spin tail is generated by recent major mergers that bias 
the distribution towards rapid rotation. Hence, the high-spin 
tail of unrelaxed haloes overlaps the distribution of quies- 
cently evolving haloes which is best fit by a log-normal dis- 
tribution. 

4 ) Alignment of halo spins and large-scale structures 

• Haloes in sheets show a strong tendency for their spin 
vector to lie in the symmetry plane of the mass distribution. 
This effect is present for all haloes but it becomes much 
more prominent for haloes with M < 10 12 /i _1 Mq. 

• For haloes in filaments, there is a slightly enhanced 
probability to find their angular momentum orthogonal to 
the filament direction, independently of mass. 

• No other significant correlation has been detected (but 
we suffer from small-number statistics in voids). 

5) Spatial correlations between halo angular momenta 

• Significant spin-spin correlations have been only de- 
tected for massive haloes in clusters. In this case, haloes 
in close pairs (separations smaller than a few Mpc) show a 
weak tendency to have antiparallel spins. 

• Alignments between spin and orbital angular momen- 
tum, however, were found to be much stronger. Regardless of 
mass and environment, spins of haloes in close pairs tend to 
be preferentially parallel to the orbital angular momentum 
of the pair. This strong effect is even enhanced for low-mass 
haloes in clusters and massive haloes in filaments. 

Our study has revealed that a number of halo proper- 
ties depend on environment. This shows that our dynamical 
classification is physical and represents a first step towards 
understanding how the galaxy formation process is influ- 
enced by large-scale structures. We will further explore the 
potential of this method in future work. 
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